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Figure  1.  Geometry  of  curve  fitting  in  cartesian  coordinates.  C  is  the  cen 
ter  of  the  gust  f ror t  feature. 

Figure  2.  Geometry  of  the  vertical  continuity  box.  Initially  determined 
size  is  dashed  and  the  solid  line  box  is  its  final  location. 

Figure  3.  Schematic  showing  detected  gust  front  and  the  two  < 'at a  windows 
over  which  we  estimate  the  horizontal  wind.  The  adjustable  pa¬ 
rameters  Ar,  A8,  and  6r  are  identified. 

Figure  4.  Illustrations  of  possible  gust  front  orientations  and  the  data 
windows  over  which  we  estimate  the  horizontal  wind. 

Figure  5.  Algorithm  output  for  the  Memphis  gust  front  of  May  28,  1985  at 

1615  CST.  Range  marks  are  at  30  km  inter  'als. 

Figure  6.  Algorithm  output  for  the  Denver  gust  front  of  July  25,  1984  at 

1732  MST  for  peak  shear  and  min  Av  thresholds  of  a)  4  m  s"*  km-* 
and  10  m  s"*,  and  b)  2  m  s"^  km"^  and  5  m  s~*.  Range  marks  are 
at  30  km  intervals.  An  "F"  indicates  a  false  feature. 

Figure  7.  Same  as  Figure  6,  except  for  the  Oklahoma  gust  front  of  April  26, 


Figure  8.  April  13,  1981  Oklahoma  gust  front  case  showing  de' acted  position 
at  a)  2111  CST  and  b)  2121  CST.  SUW  estimates  are  denoted  by  ar¬ 
rows  with  arrow  length  proportional  to  wind  speed  nd  arrow  orien¬ 
tation  indicating  wind  direction.  Mesonetwork  wine  data  are  plot¬ 
ted  as  barbs  with  flags  next  to  their  numerical  values.  Long 
flags  indicate  10  m  s“*,  short  flags  indicate  5  m  s"*.  The 
2120  CST  TTS  rawinsonde  winds  up  to  1  km  AGL  are  shown  in  (b). 

The  rawinsonde  launch  site  is  indicated  with  a  balloon.  The  line 
labeled  AA'  describes  the  location  of  the  horizontal  cross-section 
shown  in  Figure  9. 

Figure  9.  Horizontal  cross-section  along  line  AA1  shown  in  Figure  8b.  Re¬ 
flectivity  factor  Z  (dBZ),  radial  velocity  (m  s-^),  and  a  nine- 
point  average  mean  velocity  are  given.  "B“  and  "E"  denote  the 
beginning  and  ending  range  of  the  pattern  vector.  "G"  is  the 
detected  position  of  the  gust  front.  The  location  -ind  length,  in 
range,  of  each  data  window  is  shown. 

Figure  10.  Ma  9,  lc 81  Oklahoma  gust  front  case  showing  detected  position  at 
a)  0018  1ST  and  b)  0032  CST.  Symbols  are  as  defined  in  Figure  8. 

Figure  11.  April  26,  1984  Oklahoma  gust  front  case  showing  detected  position 
at  a)  2042  CST  and  b)  2052  CST.  Symbols  are  as  defined  in  Fig¬ 
ure  8.  The  1957  CST  CHK  rawinsonde  winds  up  to  1  kn  AGL  are  shown 


Figure  12.  Plot  of  the  algorithm  output  for  the  Oklahoma  gust  front  of  May  9, 
1981  along  with  the  locations  (1,2,3)  of  the  NSSL  surface  stations 
used  for  a  time  of  passage  comparison. 

Figure  13.  Same  as  Figure  12,  except  for  April  26,  1984. 

Figure  14.  Same  as  Figure  12,  except  for  April  13,  1981. 

Figure  15.  Algorithm  output  for  two  different  radars;  a)  Norman  (NRO)  and 
b)  Cimarron  (CIM),  looking  at  the  same  gust  front  (April  13, 
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1.  INTRODUCTION 

The  initial  design  and  development  of  the  gust  front  algorithm  was  done 
by  Uyeda  and  Zrnic1  (1985,  1986).  Through  limited  testing  on  several  cases, 
they  showed  that  the  algorithm  could  accurately  detect  and  track  strong  gust 
fronts  that  commonly  occur  in  Oklahoma  during  the  Spring.  However,  further 
developmental  work  and  testing  was  necessary  to  bring  the  algorithm  from  its 
semi-automatic  research  stage  to  a  routinely  functioning,  f ul ly  automated  al¬ 
gorithm  ready  for  operational  use.  The  purpose  of  this  report  is  to  document 
changes  and  additions  that  have  been  implemented  thus  far,  and  to  present  ini¬ 
tial  test  results  for  the  latest  version  of  the  algorithm.  Additional  refine¬ 
ments  will  undoubtedly  occur  after  real-time  testing  in  Denver  and  No 'man  in 
1987,  and  these  results  will  be  reported  at  a  later  date.  Readers  not  famil¬ 
iar  with  the  algorithm  are  advised  to  first  read  the  report  by  Uyeda  and 
Zrnic1 (1985). 

2.  NEW  THRESHOLDS 

The  most  important  changes  to  the  previous  thresholds  affected  those  used 
to  determine  whether  or  not  to  save  a  pattern  vector.  A  gust  front  pattern 
vector  is  a  multi -dimensional  vector  that  quantifies  radial  convergence  (a 
sequence  of  decreasing  velocities)  and  has  the  following  attributes:  azimuth, 
beginning  range  rb,  ending  range  re,  beginning  velocity  vb,  ending  velocity 
Vg ,  peak  shear,  and  location  of  the  peak  shear.  Testing  by  Uyeda  and  Zrnic' 
(1985)  indicated  that  the  best  parameter  to  use  for  determining  the  strength 
and  position  of  a  gust  front  was  the  maximum  velocity  gradient,  or  peak  shear, 
within  a  pattern  vector.  It  was  therefore  decided  to  use  the  peak  shear 
(Eq.  1  in  Uyeda  and  Zrnic',  1985)  as  the  primary  means  of  determining  whether 
or  not  to  save  a  patten  vector,  in  place  of  the  previous  high  and  low  grad- 


ient  and  "flux"  thresholds.  To  eliminate  those  vectors  having  a  high  peak 
shear,  but  low  radial  convergence  (i.e.,  a  small  velocity  difference),  a 
minimum  velocity  difference  threshold  (min  Av)  has  been  added,  with  the  velo¬ 
city  difference  given  by  Av  =  v,  -v  .  The  peak  shear  and  min  Av  thresholds 
have  initially  been  set  at  4  m  s  km_J-  and  10  m  s  ,  respectively,  which 
should  allow  for  the  detection  of  moderate  to  strong  gust  fronts  based  on 
tests  with  Oklahoma  data.  If  it  is  also  desired  to  detect  weak  gust  fronts, 
these  thresholds  will  likely  need  to  be  reduced,  noting  however,  that  as  the 
thresholds  are  reduced,  the  false  alarm  rate  will  likely  increase. 

The  criteria  for  grouping  pattern  vectors  into  possible  gust  front  fea¬ 
tures  have  also  been  changed.  Initial  tests  with  the  old  thresholds  resulted 
in  some  major  grouping  problems.  Therefore,  the  threshold  values  were  reduced 
by  about  one-half.  The  maximum  azimuthal  separation  between  v ‘Ctors  is  now 
2.2°,  so  that  a  gap  of  no  more  than  one  radial  is  allowed  if  consecutive  rad¬ 
ial  s  are  separated  by  1°.  The  maximum  range  separation  between  consecutive 
vectors  is  now  2.0  km,  with  the  range  of  a  vector  being  given  by  the  location 
of  the  peak  shear. 

In  place  of  a  height  threshold,  it  was  decided  to  process  only  two  low- 
level  scans  below  an  elevation  angle  of  1.5°  out  to  a  range  of  60  km.  Test 
results  using  up  to  four  elevation  scans  at  close  ranges  resulted  in  a  large 
increase  in  the  false  alarm  rate  for  some  cases,  with  no  significant  improve¬ 
ment  in  the  drtection  capability.  The  decision  to  process  data  to  distances 
of  60  km  was  made  to  provide  at  least  a  20  min  advanced  warning  of  gust  front 
passage  over  any  airport  close  to  the  radar  site. 

Because  it  is  distinctly  possible  that  a  feature  with  a  large  number  of 
vectors  could  in  fact  cover  only  a  small  distance  if  it  is  located  very  close 
to  the  radar,  a  minimum  length  threshold  was  added.  Initially,  a  feature  must 
have  a  minimum  of  five  vectors  to  be  saved.  Then,  after  the  length  of  a  fea¬ 
ture  is  calculated,  if  its  length  is  less  than  5  km,  it  is  discarded. 

A  list  of  the  new  threshold  values  is  given  in  Table  1. 


TABLE  l.  List  of  Thresholds. 


Threshold 
Peak  shear 

Minimum  velocity  difference 
Minimum  length 

Minim  im  number  of  vectors  in  a  feature 

Maximum  azimuthal  separation  between 
vectors 


Numerical  Valu £ 

4  m  s'1  km"-1- 
10  m  s-1 

5  k;n 
b 

2.2° 


Maximum  range  separation  between  vectors  2  km 

Minimum  reflectivity  -15  dBZ 


3.  SIGNIFICANT  ADDITIONS 

a .  Curve  Fitt ing  in  Cartesian  Coordinates 

Previously,  the  location  of  a  gust  front  was  approximated  by  fitting  a 
second-order  polynomial  r( 6)  to  the  individual  vector  locations  of  the  peak 
shears  (see  Eq .  17  in  Uyeda  and  Zrnic',  1985).  Tracking  and  forecasting  were 
performed  by  translating  the  fitted  position  according  to  the  latest  estimate 
of  propagation  speed.  While  this  procedure  performed  well  in  tracking  gust 
fronts  some  distance  away  from  the  radar,  it  may  not  be  suitable  for  fronts 
that  are  close  to  the  radar  site.  This  is  because  the  second-order  polynomial 
in  e  approximates  well  the  gust  fronts'  parabolic  shape  only  at  appreciable 
distances  f ron  the  radar.  Therefore,  it  was  decided  to  perform  curve  fitting 
and  tracking  in  cartesian  coordinates,  so  proximity  to  the  radar  would  cause 
no  problems. 

The  procedure  by  which  curve  fitting  is  done  is  shown  in  Figure  1. 

First,  a  local  coordinate  system  (x‘,  y')  is  established  with  its  origin  at 
the  center  of  a  feature.  Next,  a  line  is  formed  by  connecting  the  two  end 
points  of  the  feature,  E^  and  E£.  The  local  coordinate  system  is  then  rotated 
until  the  x'  axis  is  parallel  to  line  E-^ E 2 .  This  establishes  a  second  local 
coordinate  system  (xH,  y").  The  vector  locations  of  the  gust  front  feature 
are  then  converted  from  the  original  coordinate  system  to  the  'x",  y")  system, 
and  a  least-squares  fit  of  the  second-order  polynomial 


Figure  1.  Geometry  of  curve  fitting  in  cartesian  coordinates.  C  is  the  cen 
ter  of  the  gust  front  feature. 

b.  Vertical  Continuity 

The  init’al  version  of  the  algorithm  operated  on  only  one  elevation  scan 
of  radar  data.  This,  unfortunately,  resulted  in  a  very  high  false  alarm 
rate.  Consequently,  the  algorithm  was  changed  to  operate  on  two  low-level 
scans,  and  with  the  help  of  vertical  continuity  requirements,  to  attempt  to 
discriminate  true  features  from  false  features.  In  order  for  two  features  at 
different  elevation  scans  to  satisfy  the  vertical  continuity  requirement,  it 
is  necessary  for  the  center  of  one  feature  to  be  within  the  "vertical  contin¬ 
uity  box"  of  i ne  other  feature.  The  "vertical  continuity  box"  of  a  feature  s 
described  by  a  rectangle  extending  5  km  either  side  of  a  line  connecting  tne 


two  en  d  points  of  the  feature,  then  aligning  the  center  of  this  box  with  ere 
center  of  the  feature  (Figure  2).  The  reason  that  the  vertical  conti nuicy  oox 
extend;  the  full  length  of  a  feature  is  that,  while  the  radar  may  detect  all 
of  the  gust  front  at  one  elevation  scan,  it  is  possible  that  at  the  other 
scan,  only  a  small  oortion  may  be  detected,  due  to  either  scanning  above  most 
of  the  gust  front,  or  blockage  by  ground  clutter.  If  this  is  the  case,  and  we 
were  to  use  a  smaller  box,  it  is  also  possible  that  the  smaller  feature  would 
be  off  to  one  side  of  the  larger  feature,  and  therefore,  the  vertical  contin¬ 
uity  requirements  would  not  be  met.  Preliminary  testing  shows  that  using  two 
low-level  scans  to  identify  only  those  features  having  vertica'  continuity, 
results  in  an  ~  80%  reduction  in  the  false  alarm  rate  over  jus*  using  one  low- 
level  scan,  with  no  loss  in  detection  capability. 


Fijure  2.  7 eometrj  of  the  vertical  continuity  Doc:.  Initially  determined 

size  is  dashed  and  the  solid  line  box  is  its  fin  it  location. 


c .  Estimates  of  the  Wind  Veloci ty  Ahead  and  Behi nd  the  Gust  Pros: 

Incorporated  as  part  of  the  gust  front  algorithm  is  a  separate  algor¬ 
ithm  to  estimate  horizontal  wind,  ahead  and  behind  the  gust  front,  from  radial 
velocities  obtained  by  a  single  radar.  The  algorithm  is  a  modified  version  of 
a  procedure  known  as  the  Sectorized  Uniform  Wind  (SUW),  and  is  thoroughly  dis¬ 
cussed  in  the  report  by  Smith  (1986).  A  cursory  discussion  of  its  use  in  the 
gust  front  algorithm  is  given  in  the  following  paragraphs. 

As  implied  by  its  name,  the  SUW  algorithm  relies  on  the  assumption  of 
uniformity  of  the  wind,  over  spatial  sectors,  to  estimate  the  horizontal 
wind.  We  will  refer  to  these  sectors,  which  are  late''  to  be  specified,  as 
data  windows.  Within  these  data  windows,  the  assumed  relations  lip  between  t.he 
horizontal  wind  and  the  radial  velocity  is 

v  =  u  s  i  n  3  cos  $  +  v  cos  0  coso  f2) 

r  o  o 

where  vf  is  the  radial  velocity,  uQ  and  vQ  are  the  east-west  a-d  norm-south 

components  of  the  vector  wind,  0  is  azimuth  angle  measured  rel.tive  to  north, 

and  <j»  is  radar  elevation  angle.  For  reasons  made  apparent  by  Smith  (op. 

ci  t . ) ,  the  mean  vertical  component  of  the  wind,  wQ,  has  been  neglected  in 

(2).  The  unknown  quantities  uQ  and  vQ  are  estimated  from  a  matrix  of  redial 

velocity  measi. rements ,  dimensioned  in  range  and  azimuth,  using  the  method  of 

least-squares.  More  precisely,  the  measurements  of  radial  velocity,  v^  , 

within  a  data  window  are  regressed  on  the  spatial  varying  functions 

(sine.coso)  and  (cos0.cos$)  to  obtain  the  estimates  u  ,  v  .  It  is  assumed 
v  1  v  v  l  o’  o 

that  toe  processed  radial  winds  are  from  the  lowest  elevation  angle  of  each 
volume  scan. 

We  choose  two  data  windows  over  which  we  apply  (2)  baser  on  the  de¬ 
tected  position  of  tie  gust  front  in  range  and  azimuth  relative  to  the  radar 
(see  Figure  3).  The  detected  position  is  the  least-squares  fitted  curve 
through  the  points  of  maximum  convergent  radial  shear,  for  all  pattern  vectors 
which  describe  the  gust  front  (see  Uyeda  and  Zrnic1,  1986).  W<  see  that  these 
data  windows,  which  are  located  on  either  side  of  the  gust  fret,  are  speci¬ 
fied  by  distar ce  from  the  gust  front,  6 r,  their  range  width  Ar  (same  for  all 
radials),  and  azimuthal  width,  A9  .  All  parameters  may  be  adjustable  with 
optimal  dimensions  yet  to  be  determined.  Because  the  magnitudes  of  random 
errors  in  u^,  vQ  are  principally  controlled  by  Ao  and  inversely  related  to  it. 


a  minimum  tnreshold  on  A 9  keeps  these  estimation  errors  tolerable.  Random 
errors  in  least-squares  estimates  arise  because  errors  are  made  in  the  mea¬ 
surements  of  radial  velocity.  Typically,  we  assume  all  measurements  have  the 
same  expected  uncertainty.  Then,  the  data  window  geometry  determines  a 
scaling  factor  which  multiples  this  measurement  uncertainty  to  obtain  the  est¬ 
imate  uncertainty.  A  further  reduction  in  the  uncertainty  of  wind  estimates 
is  realized  when  radar  data  are  analyzed  over  a  range  interval  as  well  ( i . e . , 
over  Ar).  Processing  in  range  has  the  same  effect  as  averaging.  Finally,  the 
data  windows  are  displaced  in  range  by  an  amount  (6r)  from  the  detected  pos¬ 
ition  of  the  gust  front.  This  is  to  avoid  "mixing"  Doppler  velocity  data 
across  the  zone  of  maximum  radial  shear  when  applying  the  SUW  algorithm.  It 
has  been  found  that  large  departures  of  the  wind  from  uniformity  introduce 
large  systematic  errors  in  the  SUW  estimates.  An  in-depth  err >r  analysis  can 
be  found  in  Smith  (op.  cit.). 
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Fijure  3.  Schemati :  showinj  detected  just  front  and  tne  tw  >  data  windows 
over  which  we  estimate  the  horizontal  wind.  The  adjustable  pa¬ 
rameters  Ar,  A9,  and  <5r  are  identij  ied. 


The  designated  locations  of  th  ;se  data  windows  relative  to  the  detected 
gust  front  can  best  be  illustrated  wi  :h  examples.  Figure  4  shows  schematics 
illustrating  possible  orientations  of  gust  fronts  relative  to  the  radar  (lo¬ 
cated  in  the  center  of  the  figure)  an  i  the  data  windows  where  Doppler  data  are 
processed  to  obtain  mean  wind  estimates.  (The  present  version  of  the  gust 
front  algorithm  does  not  have  the  capibility  of  measuring  azimuthal  shear  am 
hence  detecting  gust  fronts  parallel  :o  radar  radials.  Consequently,  we  will 
only  concern  ourselves  with  those  gust  fronts  with  orientations  nearly  perpen¬ 
dicular  to  radar  radials.)  For  example,  in  Figure  4a,  the  center  azimuth  of 
the  gust  front  is  at  a  range  nearer  the  radar  than  its  end  points.  In  addi¬ 
tion,  let's  assume  that  the  azimuthal  width  is  found  to  be  less  than  the  mini¬ 
mum  threshold  we  imposed  on  A 9.  Themfore,  data  at  an  equal  number  of 


(c)  I  (d) 


Fijure  4.  Illustrations  of  posy.  t>le  just  front  orientations  and  the  data 
windows  ooer  which  estimate  the  horizontal  wind. 


a/.iinutus  or-,  either  side  of  its  end  points  are  used  to  construct  an  appropriate 
dita  w  ndow  (i.e.,  so  that  the  azimuthal  extent  of  the  data  window  is  at 
least  AQwide).  In  Figure  4c,  the  azimuthal  extent  of  the  gust  front  is  below 
tie  imposed  minimum  threshold  and  its  orientation  is  such  that  the  range  at 
one  end  point  is  closer  to  the  radar  than  its  center  range  and  other 
endpoint.  In  this  case,  the  data  window  is  extended  in  azimuth,  at  constant 
range,  fron  the  end  point  of  the  gust  front  closest  to  the  radar.  Fur  he^ 
examples  (Figures  4b  and  4d)  show  differing  orientations  of  gust  front  and 
the  resulting  data  windows.  The  important  point  is  that  when  the  gust  rronts 
are  of  azimuthal  width  at  least  equal  to  the  minimum  threshold  (i.e.,  ia), 
data  over  the  entire  az  muthal  width  af  the  gust  front  are  used  to  estimate 
the  mean  horizontal  winds;  otherwise,  the  azimuthal  width  of  the  analysis 
region  must  be  extended. 

4.  PRELIMINARY  TEST  RESULTS 

a.  Test  Results  for  the  Gust  > ront  41 gori tnm 

To  he  classified  as  a  si gni f i <  amt  gust  front,  a  radial  velocity  differ- 
ence,  Av,  of  at  least  S  m  s'1  and  a  radial  shear  of  at  least  2  m  s'1  km-1  must 
extend  over  an  azimuthal  length  of  5  km  or  more,  and  persist  at  least 
10  minutes.  We  further  classify  the  strength  of  a  gust  front  based  on  the 
average  radial  velocity  difference  o-  ar  its  azimuthal  extent  as  follows: 
weak,  Av  =  b-9  n  s-^;  moderate,  Av  =  10-14  m  s"*;  strong,  Av  =  15-25  in  s~^; 
and  severe,  Av  >  25  m  s"*.  The  dete-  nination  of  correct  and  false  algorithm 
detections  wa^  based  on  the  single  D-ppler  radar  reflectivity  and  velocity 
fields.  Most  false  detections  occur  ad  in  the  groan1  clutter  as  the  ’ow-leve" 
winds  and  grr  nd  clutter  interacted  )  produce  false  shear  zones.  Occasion¬ 
ally,  a  non-gust  front  shear  zone  wa  detected  in  an  area  of  active  convectio- 
(e.g.,  second. ry  surge  or  obstacle  f.jw).  These  were  considered  as  neutral 
detections,  aid  not  counted  as  eithe  correct  or  fals?  detections.  Because 
the  purpose  of  the  algorithm  is  to  d  tect  and  track  only  those  gust  fronts 
that  are  aoprr, aching  or  receding  fro  the  radar,  gust  fronts  with  orientation 
nearly  parallel  to  the  radar  beam  wi  1  likely  not  be  detected,  or  may  be 
classified  as  weaker  than  they  actua  ly  are. 

Preliminary  tests  were  run  fc  the  latest  version  of  the  gust  front 
algorithm  on  t.he  following  cases  fro  the  Nationa,  Severe  Storms  Laboratory 
(NSSL),  FAA/Li ncoln  Lab,  and  the  Nat  onal  Center  For  Atmospheric  Research 


(NCAR).  They  are: 


4/13/81,  2105-2127  CST,  five  volume  scans 

1)  3  NSSL  cases  5/9/81,  0006-0048  CST,  four  volume  scans 
4/26/84,  2028-2052  CST,  five  volume  scans 

Data  for  each  of  these  cases  encompassed  strong  to  severe  gust  fronts 
(maximum  radial  velocities  behind  the  gust  front  were  25-30  m  s"*)  whose 
lengths  varied  from  15  to  55  km.  Using  the  thresholds  shown  in  Table  1,  the 
algorithm  worked  very  well  in  detecting  and  tracking  the  gust  fronts.  For  the 
4/13/81  and  5/9/81  cases,  there  were  no  false  alarms.  For  the  4/26/84  case, 
there  was  either  zero  or  one  false  alarm  per  volume  scan,  with  the  average 
being  0.6  per  volume  scan.  The  algorithm  outputs  for  these  three  cases  were 
in  agreement  with  the  results  of  Uyeda  and  Zrnic'  (1985). 

2)  1  FAA/Lincoln  Lab  Operational  Weather  Study  (FLOWS)  case: 

5/t8/85,  1558-1651  CST,  eleven  volume  scans 

This  case  started  out  as  a  small,  narrow,  weak  gust  front  (maximum  rad¬ 
ial  velocities  behind  the  gust  front  were  6-9  m  s"*).  It  later  developed  into 
a  broader,  more  diffuse  outflow.  In  order  to  detect  this  feature,  it  was  nec¬ 
essary  to  reduce  the  minimum  velocity  difference  threshold  (min  Av)  to  5  m 
s"1,  the  peak  shear  threshold  to  2  m  s“*  km’*,  and  the  minimum  length 
threshold  to  3  km.  Even  then,  detection  occurred  only  about  half  of  the  time, 
due  to  the  weak  nature  of  this  front.  Surprisingly,  even  with  this 
significant  reduction  in  the  thresholds,  there  were  no  false  alarms.  An 
example  of  algorithm  output  for  one  of  the  times  when  the  gust  front  was 
strongest  is  shown  in  Figure  5. 

3)  1  NCAR  Classify,  Locate,  and  Avoid  Wind  Shear  (CLAWS)  case: 

7/25/84,  1621-1728  MST,  fourteen  volume  scans 

This  case  consisted  of  one  long,  weak  to  moderately  strong  convergence 
line  and  several  small,  weak  to  moderately  strong  gust  fronts  (maximum  radial 
velocities  behind  the  gust  front  of  6-12  m  s"*).  Initial  tests  with  the 
thresholds  in  Table  1  (set  for  the  strong  Oklahoma  cases)  resulted  in  the  de¬ 
tection  of  about  half  of  the  convergence  line  during  the  times  it  attained  a 
moderate  strength;  weaker  portions  were  not  detected.  None  of  the  small  scale 
gust  fronts  were  detected,  and  there  were  no  false  alarms.  The  thresholds  for 
min  Av  and  peak  shear  were  then  reduced  to  5  m  s’*  and  2  m  s’*  km’*,  respec¬ 
tively,  duplicating  two  of  the  reduced  thresholds  in  the  FLOWS  case.  With 


this  change  we  were  able  to  detect  all  of  the  converg *nce  line  ;nd  the  small 
scale  gust  fronts.  However,  this  brought  about  an  increase  in  the  false  alarm 
rate  to  an  average  of  one  per  volume  scan.  Figure  fi  shows  algorithm  output 
for  both  sets  >f  thresholds  near  the  time  of  maximum  shear  alone  the  conver¬ 


gence  line. 


28  MAY  1985  1615  CST 


Figure  5.  Algorithm  output  for  the 
Memphis  gust  front  of  Ma /, 
28,  1985  at  1615  CST. 
Range  marks  are  at  20  km 
intervals. 
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Algorithm  output  for  t h-  Denver  gust  front  of  Jut.j  , 
1702  M^T  for  peak  shear  ind  min  A V  thresholds  of  a1  ; 
and  10  m  s~  ^ ,  and  b)  2  n  s~  ^  km~  ^  an!  5  m  1 ,  Range, 
at  50  m  intervals .  An  ’F"  indicates  a  'alse  feilur-t. 
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Since  t.hc  nhovf  results  from  the  FLOWS  and  NCArt  cases  showed  that  a  re  - 


duct  i  oi  in  the  thresholds,  in  order  to  detect  smaller,  weaker  “eatures,  did 
not  result  in  a  large  increase  in  false  alarms,  a  check  was  nade  to  see  the 
effect  such  a  reduction  would  have  on  one  of  the  stronger  Oklahoma  cases. 
Therefore,  a  second  test  using  the  4/26/34  case  was  run  with  the  nin  pv 
threshold  at  -  ~  s""  and  the  peak  shear  threshold  at  2  m  s-^  km'".  Tnis  v~e, 
however,  a  larje  increase  in  false  al  trms  did  result,  with  the  average 
increasing  from  0.6  per  volume  scan  to  3.6  per  volume  scan.  What  tnis  seems 
to  show  is  that,  for  the  dynamically  stronger  weather  environments  associated 
with  long,  intense  gust  fronts,  it  i'  necessary  to  retain  the  thresholds  at 
their  iriginal  values  in  Table  1  to  keep  the  false  alarm  rate  low,  whereas, 
for  the  dynamically  weaker  weather  ewironments  associated  with  small-- n,  less 
intense  gust  fronts,  it  may  be  possi*  le  to  reduce  the  thresholds  and  •  .ill 
have  a  low  false  alarm  rate.  Figure  7  shows  algorithm  output  for  both  sets  of 
thresholds  at  203/  CSf. 
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What  tie  above  results  show  is  that  the  performance  of  tne  ij-.r-t  '  ront. 
algori thm  depends  largely  on  what  type  of  feature  it  is  to  detect,  which  in 
turn  determines  the  values  of  the  thresholds  to  oe  us»d.  With  the  orijina! 
thresholds,  the  algorithm  works  very  well  at  detecting  the  larger  scale, 
strom  gust  fronts  associated  with  fast-moving  sgu.al  1  lines  and  bow  eihoes; 
false  alarms  are  few,  if  any.  As  the  thresholds  are  reduced  to  detect  weak  to 
moderate  gist  fronts,  the  number  of  false  alarms  generally  goes  up.  There¬ 
fore,  1  decision  must  be  made  as  to  woat  the  proper  balance  should  be  between 
trying  to  detect  increasingly  weaker  just  fronts  with  a  correspondi ng  poten¬ 
tial  increase  in  the  false  alarm  rate,  bearing  in  mind  that  lowering  the 
thresholds  in  the  more  dynamic  weather  situations  is  much  more  deleterious  to 
the  false  alarn  rate  than  similar  reductions  of  the  thresholds  in  less  dynamic 
weather  conditions. 

b •  Test  Results  for  the  Sectorl zed  Uniform  Wind  A1 gori thm 

Three  gust  front  cases  in  central  Oklahoma  were  chosen  to  test  tne  S U 
algorithm;  April  13,  1981,  May  9,  19  l,  and  April  26,  1984.  These,  as  well  as 
some  other  cases,  are  documented  in  .‘tail  in  the  report  by  Uyeda  and  Zrnic' 
(1985).  Doppler  velocity  displays  f  ~  several  of  the  analysis  times  can  be 
found  in  the  above  cited  report. 

We  compared  the  wind  output  f-  >m  the  SIJW  algorithm  with  other  supple¬ 
mental  wind  data  such  as  that  availa  le  from  Stationary  Automated  Mesonetworr 
(SAM)  sites.  When  available,  winds  rom  rawinsonde  and  those  measured  by  an 
instrumented  tower  were  used  to  jtidg-  SlJW  results.  Comparisons  of  rawinsonde 
data  with  SIJW  estimates  are  best  mad'  ahead  of  the  gust  front  since  both  these 
measurements  represent  to  some  exten  the  ambient  environmental  wind.  Fur¬ 
thermore,  the  wind  ahead  of  the  fron'  tends  to  be  fairly  uniform  and  thus  com 
plies  with  thf  basic  assumption  of  t  e  SOW  algorithm. 

For  thf se  three  analyzed  case  ,  the  following  values  for  the  parameter 

descriiing  t h.  data  window  size  and  ocation  were  used:  Ad  -  40°,  Ar  ranged 
betwee  i  2.25  in  and  3.15  km,  and  -5r  anged  between  459  m  and  630  m,  de;  ending 

on  the  radar  setup.  The  radial  spac  ng  of  the  velocity  data  was  either  150  m 

corresponding  to  the  first  set  of  va  ues  for  at  and  or,  or  210  m  corresponiin 
to  the  second  set  of  values. 
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Apri  1  13,  1931 : 


Figures  3a  and  8b  illustrate  the  wind  comparisons  made  at  2111  CsT  and 
’123  (1ST,  respectively.  At  approximately  2120  CST,  a  rawinsomie  was  released 
r  rom  lut-tle,  Oklahoma  (118).  IIS  is  located  almost  due  west  of  the  NV>i  rmla'- 
^NKO'  near  Mi  vin  range.  The  rawinsondo  winds  up  to  1  km  above  ground  level 
(Alii.)  are  show  i  on  Figure  8b.  Mesonetwork  wind  data  recorded  at  analysis  tine 
are  superimposed  on  the  figures. 


At  2111  CST,  general  agreement  between  SUW  estimated  wind  direction  and 
those  of  the  mesonetwork  data  is  apparent  (Figure  8a).  On  the  other  hand,  the 
wind  soeeds  for  the  SUW  estimates  ter  i  to  be  higher.  This  is  most  evi  ient 
ahead  (east)  of  the  gust  front.  However,  the  radar  data  used  to  compute  the 
SUW  estimates  ire  from  near  750  m  AO!  ,  whereas  the  mesonetwork  data  are  reore- 
sentative  of  winds  near  the  ground.  Furthermore,  examination  of  Doppler  dis¬ 
plays  indicate  radial  velocities  nea  11  m  s-^  immediately  ahead  of  the  gust 
front  and  near  19  m  s“^  immediately  *hind.  Thus  the  SUW  estimates  of  wind 
magnitude  agree  reasonably  well  with  the  largest  measured  radial  winds. 


At  2123  CST  (Figure  8b),  consid  rable  differences  are  observed  not  only 
between  SUW  estimated  winds  and  the  AM  winds  but  also  between  the  SAM  winds 
and  those  measured  from  the  rawinsonde.  For  example,  the  surface  measurement 
at  the  rawinsonde  site  differs  nearly  10  m  s_1  in  wind  speed  from  a  SAM  mea¬ 
surement,  even  though  this  network  station  is  located  only  a  few  kilometers  to 
its  north.  In  general,  however,  the  network  stations  consistently  indicate 
weak,  westerly  winds  in  advance  of  V  ?  gust  front,  which  veer  to  the  northwest 
behind  the  gust  front  with  speeds  ne  r  15  m  s"^. 


'  rom  inspection  of  the  Doppler  elocity  display  at  this  same  time  'see 
Figur  13  in  I yeda  and  Zrnic',  1985)  we  conclude  that  the  SUW  estimate  of 
wind  rirectior  ahead  of  the  front  is  reasonable.  It  not  only  agrees  well  wit' 
the  r-.winsondf  winds  at  the  same  alt  tude  (radar  beam  height  about  400  m  ACL), 
but  the  line  .  f  ^ero  Doppler  velocit  is  aligned  nearly  perpendicular  to  the 
estimated  win-  direction.  From  the  agnitudes  of  the  Doppler  winds  ahead  of 
the  front,  we  can  not  easily  assess  he  accuracy  of  the  estimated  wind  magni¬ 
tude  but  we  note  general  consistency  between  the  two;  the  Doppler  winds  are 
relatively  weak  (magnitudes  less  tha  10  m  s~M  and  so  is  the  SUW  wind.  Be¬ 
hind  the  front,  mean  wind  direction  s  difficult  to  determine  because  of  the 
diver-iing  nati  re  of  the  outflow.  Ho  ever,  velocity  aliasing  nas  occurred,  in- 
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7igur-.  8.  April  13,  1.981  Oklahoma  gust  front  case  shoving  detected  position 
at-  a)  2111  CST  and  b )  2123  CST.  SUW  estimates  are  denote 1  bg  ar¬ 
rows  with  arrow  length  proportional  to  wind  speed  and  arrow  orien¬ 
tation  indicating  wind  direction.  M isonetwork  wind  data  are  plot¬ 
ted  as  barbs  with  flags  next  to  their  numerical  values.  Long 
flags  indicate  10  m  s~ 7  short  flags  indicate  5  m  s~ ^ .  The 
2120  CST  TTS  rawinsonde  winds  up  to  1  km  ACL  are  shown  in  (by. 

Tl  e  rawinsonde  launch  site  is  indicated  with  a  balloon.  The  line 
V  baled  AA'  describes  the  location  of  tne  horizontal  cross-section 
3'  O'jn  in  Figure  9. 


'* 

'* 


Figure  9  illustrates  a  horizontal  cross-section  through  the  gust  front  nr, 
2L23  CST  showing  reflectivity  factor  and  radial  velocity  along  the  310° 
radial.  Radar  beam  elevation  angle  is  1.2°.  The  detected  position  of  the 
gust  front  ( <3 ) ,  and  the  beginning  (3)  and  ending  '?.)  °ange  of  the  pattern 
vector  are  indicated.  The  range  limits  of  the  two  data  windows  are  also 
shown.  As  i I  Instated,  the  region  outlined  as  the  pattern  vector  is  cnaracter- 
izeii  b /  convergent  radial  shear. 


Figure 


horizontal  cross-section  a  rng  line  AA’  shown  in  Figure  So.  Re¬ 
flectivity  factor'  Z  (d3Z),  r adial  velocity  (m  s~^)}  and  a  ninc- 
ooint  averaje  mean  veiocit.  are  given.  "B"  and  "S"  denote  sac 
be j  nniny  and  ending  range  )f  the  pattern  vector .  "G"  is  tne 
let  cted  position  of  the  j  tt  front.  Tne  location  and  lengtn ,  In 
ran  .e,  of  each  <iata  window  is  shown. 
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An  observer;  characteristic  of  many  tnunderstorm  outflows  is  a  w’nd  max' - 
num  immediately  oehird  the  gust  front  [Figures  12  and  13  in  Eilts  '1936a)  are 
good  examples].  Goff  (1976)  attributes  this  to  a  solenoidal  circulation  near 
the  boundary  b  :tween  cold  outflow  and  warm  inflow,  which  accelerates  the  flow 
in  the  cold  ai  ~  immediately  behind  the  gust  front.  In  Figure  9,  tneso  pear, 
inbound  winds  occur  where  the  radial  shear  is  zero,  some  6-7  km  behind  the 
gust  front.  Behind  this  wind  maximum,  there  is  a  region  of  moderate  di verged 
radial  shear,  beyond  which  the  winds  remain  fairly  steady  but  still  relatively 
strong.  As  illustrated,  the  data  window  behind  the  gust  front  encompasses, 
for  the  most  part,  this  region  of  maximum  winds.  Therefore,  it  is  not  sur¬ 
prising  that  the  SA  1  winds,  measured  well  behind  the  gust  front,  are  weaker  in 
magnitude.  Furthermore,  as  pointed  out  earlier,  the  SUW  estimate  represents 
the  winds  near  400  m  AGL.  Eilts  (1986b)  observed  that  measured  outflow  winds 
at  heights  a bo /e  200  m  tend  to  be,  on  the  average,  1.6  times  larger  than  winds 
measured  at  the  surface.  If  we  multiply  the  surface  wind  speeds  by  a  factor 
of  1.6,  they  agree  quite  well  with  the  speeds  measured  aloft. 

May  9,  1981: 

figures  11a  and  10b  illustrate  the  SIJW  output  for  0018  and  0032  CST,  re¬ 
spectively.  Very  little  supplemental  wind  information  is  available,  primarily 
due  to  the  far  range  of  the  gust  front  (i.e.,  it  is  outside  tne  range  of  toe 
surface  mesonetwork ) .  In  addition,  there  was  no  proximity  rawinsonde  re¬ 
lease.  Therefore,  much  of  our  assess  tent  of  the  accuracy  of  the  S:JW  estimates 
relies  on  inferences  drawn  from  the  imppler  velocity  displays. 

At  the  earlier  time  (0018  CST),  >IJW  output  indicates  a  pronounced  wind 
shift  across  the  front  of  nearly  90°  '196°/13  in  s“^  to  282°/17  m  s“*).  Unfor¬ 
tunately,  all  available  surface  observations  occur  well  to  the  east.  We  did 
find,  however,  that  the  line  of  zero  loppler  velocity  ahead  of  the  front 
agrees  quite  well  with  the  SIJW  estimc  .e  of  wind  iirection.  Fortunately,  by 
0032  CST,  the  gust  front  position  was  near  the  western  edge  of  the  SAM  net¬ 
work.  At  this  later  time,  the  SUW  estimate  ahead  of  the  gust  front  co'  pares 
extremely  well  with  the  surface  obser  /ations  both  in  speed  and  direction. 

There  is  no  surface  wind  information  ivailable  for  comparison  behind  the  gust 
front.  The  0033  CST  Doppler  velocit,  display  (see  Figure  21  in  'Jyeda  and 
Zrnic1,  1985)  indicates  near  20  m  s'  radial  winds  immediately  behind  the  gus  + 


1/ 


front.  The  strength  of  these  winds  <ire  in  close  agreement  wifi  the  hmi.mntal 
wind  i-itimated  from  the  radial  velocity  measurements. 
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9  MAY  1981 


~80  km 


Figure  10.  May  9,  1981  Oklahoma  gust  front  case  showing  detected  position  a: 

a)  0018  GST  and  b)  00'6v.  CST.  Symbols  are  as  defined  in  Figure  8. 


April  26,  1934: 

Like  the  April  13,  1981,  gust  front  case,  there  is  an  abundance  of  sj?- 
plemental  wind  data  on  this  day.  Fi.ures  11a  and  lib  illustrate  the  wind  com¬ 
parisons  for  042  and  2062  CST,  resp  ctively.  A  rawinsonde  release  from  CHK 
(Chickasha,  Oklahoma)  occurred  at  20  7  CST.  The  CHK  launch  site  is  located 
approximately  40  km  southwest  of  Nor  ;an  (NRO).  The  rawinsonde  wind  data  is 
plotted  in  Figure  lib  for  the  lowest  kilometer. 
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April  26 ,  1984  Oklahoma  ust  front  ease  shoeing  detected  position 
at  a)  2042  CST  and  b)  20-2  CST.  Symbols  are  as  defined  in  Fig¬ 
ure  8.  The  1967  CST  CHK  rawinsonde  uinds  up  to  1  km  ACL  are  shoun 
in  (b). 
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At  20^2  CST,  the  leading  e  Ige  of  the  gust  front  is  located  around  4: 
northwest  of  NRO,  well  within  the  SAM  network.  Wo  see  excellent  agreement  be¬ 
tween  the  SIJW  estimated  wind  directions  and  those  measured  at  the  surface. 
Again  the  magnitudes  of  the  SUW  estimates  are  slightly  larger.  The  line  of 
zero  Doppler  velocity  on  the  203/  CST  display  indicates  south-southwest  winds 
ahead  of  the  gust  front,  which  is  in  excellent  agreement  with  the  SUW  estimate 
of  direction. 

Comparisons  made  at  2052  CST  (Figure  lib)  show  generally  small  differ¬ 
ences.  The  SUW  estimate  ahead  of  the  front  indicates  a  relatively  strong, 
southerly  wind  (20  m  s-^  at  179°).  '  oe  rawinsonde  measured  winds  near  the 

same  height  (approximately  300  m  ACL'  are  more  southeasterly  and  much  weaker 
(10  m  less),  which  is  expected  because  the  SUW  estimates  are  derived  f-om 
radar  lata  very  near  to  the  front  wh.  re  accelerations  of  the  environmental 
flow  recur.  Rehind  the  gust  front,  me  SUW  estimate  indicates  a  westerly  wind 
of  13  :i  s“*,  whereas  the  surface  val-.es  are  more  northwesterly  and  quite  a  bit 
weaker . 

In  addition  to  the  rawinsonde,  there  were  wind  measurements  from  the  in¬ 
strumented  KTVf-TV  tower  located  in  Oklahoma  City,  some  36  km  north  o* 

Norman.  The  t  >wer  is  instrumented  at  seven  heights,  the  highest  level  is  at 
444  n  AGL.  On.  minute  averages  of  wind  speed  and  direction  at  the  seven  re¬ 
cording  levels  are  shown  in  Table  2  for  2042  CST  (several  minutes  prior  to 
gust  front  passage)  and  2100  CST  (ap>  roximately  4  minutes  after  the  gust  front 
passed  the  tower). 

The  wind  speeds  recorded  at  the  tower,  above  177  m  and  prior  to  gust 
front  passage,  are  in  oetter  agreement  with  the  prefrontal  wind  speeds  esti¬ 
mated  from  the  Doppler  rudial  winds  ' han  those  measured  by  the  rawinsonde  near 
the  same  heights  (see  Fig.  11a  and  b  .  Furthermore,  the  wind  directions  mate, 
the  SUW  estimates  extemely  well.  Fo  lowing  gust  front  passage  of  the  tower, 
around  2100  CST,  an  overall  westerly  wind  is  indicated,  which  is  in  better 
agreement  with  the  SUW  estimates  of  ,'ind  direction  than  those  from  SAM  mea¬ 
surements.  This  is  not  surprising  s  nee  the  tower  recorded  winds  at  225  in  and 
444  m  are  near  the  heights  represent  d  by  the  radar  data  used  in  tne  S'JW  anal¬ 
ysis.  What  is  somewhat  surprising  i  that  the  SUW  estimates  of  wind  speed 
overestimate,  by  50%,  the  wind  speed  measured  at  the  upper  levels  of  the 
tower.  This  might  well  be  explained  by  the  fact  that  the  tower  is  located  a 
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minimi  i  of  M)  km  from  Che  centroid  of  the  gust  front  where  the  SUW  estimates 
are  mo  it  representati ve.  We  believe  that  the  SUW  estimates  are  reasonable 
because  the  Ooppl^r  velocity  fields,  at  the  times  of  the  SUW  analyses,  indi¬ 
cate  winds  behind  the  gust  front  much  stronger  than  the  tower  winds  recorded 
at  2100  CSr. 


fable  9 .  ’owe  •  data  at  2042  and  2100  CST  for  the  Apri 1  26,  1904 
Oust  Front  Case. 


2042  csr 


2100  CST 


Lev-  1 


17/  m 
266  m 
444  m 


direction  (3)/Speed  (m  s'1) 


168/13.9 
1/0/lb. 6 
175/15.9 
1/1/15.9 
169/19.3 
171/23.7 
170/27.0 


direction  (°)/Speed  (m  s 


274/11.1 

278/11.3 

278/12.0 

276/12.3 

274/13.4 

272/12.8 

244/10.0 


c.  Comparison _of  SOW  Estimate  With  the  Propagation  Speed  of  the  ujst 
Front 

The  propagation  speed  of  the  ust  front  on  readily  be  determined  if 
its  position  is  known  at  two  consecutive  times,  tp  and  t£»  Because  parts  of 
the  gust  front  may  move  relative  to  --ach  other,  thus  causing  a  change  in 
shape,  a  mean  speed  is  calculated  ov  r  some  time  interval  At  =  1 2~t ^  from  the 
displacement  of  its  centroid.  The  o  ntroid  is  defined  by  the  gust  front's 
center  range  and  azimuth  (Eqs.  5  and  5  in  Uyeda  and  Zrnic',  1985). 

The  three  analyzed  cases  present  opportunities  to  compare  the  propaga¬ 
tion  speed  with  the  SUW  estimates  of  wind.  Various  relationships  between  the 
propagation  speed  of  the  gust  front  and  the  mean  wind  components  ahead  and  be 
hind  the  gust  front,  and  normal  to  tie  frontal  boundary,  have  been  proposed. 
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Goff  (19/6)  suggests  that  the  winds  in  the  warm  air  (i.e.,  ah  >od  ot  the  gu  ,t 
front)  have  little  influence  and  that  the  propagation  speed  i  m  he  described 
>piite  accurately  as  two-thirds  of  the  peak  wind  measured  behind  the  gust 
front.  In  Table  3,  we  list  the  estimated  propagation  speed  along  with  the  SIM 
estimates  of  wind  behind  the  front. 

Table  3.  Gust  Front  Propagation  Speeds  Obtained  From  Two  Consecutive  (In 

”ime)  Centroid  Positions  Along  With  SUW  Estimates  of  the  Wind  Beninr 
the  Gust  Front  for  the  Earlier  of  the  Two  Times. 


HATE 

CENTROID 

TRACKING 

SUW  ESTIMATE 

I3EHI  ID  GJST  FRONT 

DIRECTION  (°) 

SPEED  (m 

s'1)  DIRECTION 

(°)  SPEED  (m  s'1; 

4/13/31 

288 

18 

289 

18 

b/(19  31 

303 

19 

232 

17 

4/26/84 

28 

12 

284 

16 

For  the  two  cases  in  1981,  the  SUW  estimated  winds  and  the  propagation 
speeds  agree  reasonably  well,  i.e.,  no  scaling  factors  are  required  to  match 
the  SUW  estimate  with  the  speeds  determined  from  centroid  tracking.  A  large 
discrepancy  exists  for  the  1984  case,  however.  Although  the  displacement  of 
the  centroid  indicates  a  movement  to  the  south-southwest.  Figures  3a  and  8o 
clearl/  show  movement  to  the  southeast.  We  believe,  in  this  case,  that  the 
SUW  estimate  of  wind  agrees  better  with  the  actual  propagation  speed  than  that 
obtained  f’-om  centroid  tracking.  A  problem  that  plagues  the  tracking  of 
larger  gust  fronts,  such  as  the  one  rn  4/26/84,  is  that  the  centroid  position 
can  change  dramatically  in  the  direction  along  the  front,  for  relatively  small 
changes  in  gust  front  position.  Furthermore,  if  at  some  time  part  of  the  gust 
front  becomes  aligned  along  radar  racial s,  propagation  speeds  obtained  from 
tracking  centroid  positions  may  be  totally  unrepresentative  of  the  true  propa¬ 
gation  speed.  More  reliable  estimates  of  propagation  speed  might  be  obtained 
if  relationships  such  as  those  found  by  Goff  (1976)  are  used. 
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fl.  G^erul  Comments  on  the  SUW  Algoritnm 

One  nf  the  objectives  of  incorporating  the  S'JW  algoritnm  into  the  gust 
front  algorithm  is  to  estimate  the  wind  shear  across  the  front.  As  is  well 
known,  this  is  invaluable  information  to  pilots  because  it  guantiMes  the 
hazard  and  thus  aids  in  aircraft  safety.  Surface-based  sensors  can  always 
provide  wind  information  near  airport  runways.  However,  Doppler  radar  can 
supply  the  winds  some  distance  from  the  airport,  possibly  well  in  advance  of 
gust  f~ont  passage.  This  wind  information  can  aid  air  traffic  controllers  i- 
decisi  in  makinj  for  changing  runways.  With  the  three  cases  presented,  a 
sample  of  the  performance  of  the  SUW  algorithm  was  given.  Quantitative  com¬ 
parisons  were  limited  by  the  availability  of  wind  data  from  the  SAM  network, 
rawinsondes,  and  the  tower. 

It  is  important  to  consider  situations  which,  when  they  arise,  may  make 
it  difficult  to  estimate  the  horizontal  winds  from  single  Doppler  radar  mea¬ 
surements.  As  is  well  known,  the  minimum  height  at  which  the  gust  front  can 
be  detected  by  radar  depends  on  several  factors,  many  of  which  the  radar  oper¬ 
ator  has  little  or  no  control  over.  For  instance,  the  earth's  curvature  sets 
a  lower  limit  to  the  observation  height  at  any  range.  At  near  ranges,  the 
earth's  curvature  may  not  be  of  great  concern.  However,  the  gust  front  can 
become  “lost"  in  the  ground  clutter.  Ground  clutter  cancellors  should  allevi¬ 
ate  this  problem  somewhat  and  allow  tracking  very  near  the  radar  site.  It  is 
at  these  close  ranges  that  the  accuracy  of  SUW  estimates  degenerates  i£  t^ac*- 
ing  is  accomplished  by  increasing  the  elevation  angles  of  observation.  S~mt 
(op.  cit.)  shows  that  random  errors  in  wind  estimates  are  inversely  reated 
to  cos's*.  In  addition,  artifacts  such  as  range  overlaid  echoes  and  velocity 
aliasing  need  to  be  considered.  Velocity  aliasing  may  make  SUW  estimates  use- 
1  ess . 

Furthermore,  we  need  to  consider  the  parameters  which  describe  the  dc.ta 
window  and  its  position  relative  to  the  detected  gust  front.  Although  we  used 
40°  a.  the  minimum  threshold  for  AO,  we  can  tolerate  azimuthal  widths  as  small 
as  30  .  Random  errors  in  SUW  estimates  do  not  increase  substantially  for  this 
smaller  value.  One  advantage  of  processing  Doppler  data  over  smaller  sectors 
is  that  more  reliable  estimates  of  the  winds  behind  the  smaller  gust  fronts 
may  be  possible,  particularly  for  those  gust  fronts  which  ar>  of  azimuthal 
width  less  than  the  minimum  threshold  imposed  on  A3.  This  i  because  less 
data  "outside1  the  bounds  of  the  gust  front  would  be  used.  In  addition,  the 


effects  of  oorniri i  form t i os  of  the  wind  on  SUW  estimates  tend  to  he  le->s  ov  »r 
these  smaller  sectors  'see  Smith,  op.  cit.). 

Although  the  range  width  Ar  was  chosen  arbitrarily,  we  nee  1  to  place  an 
j  yy-t1'  ’  i  m  i  r  on  this  value  as  well.  As  the  elevation  angle  increases,  the 
-ada'-  .earn  traverses  an  increasing  increment  in  height  for  a  given  range  in- 
t  -vu'.  Conversely,  f io r  a  given  elevation  angle,  increasing  tne  "ange  inter¬ 
val  increases  the  height  increment.  This  may  present  problems  if  vertical 
wind  shear  is  large.  Vertical  shear,  as  well  as  other  nonuniformities  such  as 
secondary  surges  and  downdrafts,  contaminate  SUW  estimates.  Therefore,  the 
range  widths  must  not  he  excessive;  otherwise,  the  basic  assumption  of  a  uni¬ 
form  wind  over  the  analysis  area  are  likely  to  be  violated.  We  believe  range 
widths  less  than  10  km  are  acceptable. 

Finally,  the  range  displacement  of  the  data  windows  from  the  detected 
range  of  the  gust  front  6r  may  need  to  be  adjusted.  As  discussed  earlier,  the 
strongest  winds  are  displaced  somewhat  behind  the  gust  front.  Presently,  the 
data  window  behind  the  gust  front  encompasses,  for  the  most  part,  the  strong¬ 
est  winds.  It  may  be  necessary  to  use  two  data  windows  so  that  we  not  only 
estimate  the  strongest  winds,  but  also  the  more  sustained  winds  well  within 
the  thunderstorm  outflow.  The  data  window  ahead  of  the  gust  front  may  need  to 
be  further  displaced  so  that  it  is  situated  some  distance  away  from  the  up¬ 
draft  region,  which  generally  occurs  in  the  warm  air  as  it  is  lifted  over  the 
cold  outflow.  Goff  (op.  cit.)  has  found  this  updraft  region  to  be  narrowest 
in  the  horizontal  direction  near  the  surface,  typically  no  more  tnan  a  kilome¬ 
ter  wide;  the  updraft  region  then  broadens  with  height  to,  at  most,  several 
kilometers.  Displacement  of  the  data  window  by  at  least  5  km  would  assure 
tnaf  foe  estimated  winds  are  not  influenced  by  the  pre-frontal  updraft. 

5.  VERIFICATION  OF  GUST  FRONT  ALGORITHM  OUTPUT 

a.  Comparison  of  the  Gust  Front  Position  with  Surface  Stations 

To  determine  how  well  the  location  of  the  gust  front,  as  indicated  by 
the  algorithm,  compared  with  its  location  as  indicated  by  surface  stations, 
algorithm  output  for  the  three  NSSL  cases  was  compared  with  the  surface  wind 
direction  obtained  from  several  NSSL  surface  stations.  Since  it  is  not  possi¬ 
ble  to  determine  the  exact  location  of  a  gust  front  from  a  group  of  scattered 
surface  stations,  we  decided  to  compare  the  time  of  passage  of  the  wind-shift 
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line  over  the  available  surface  sites.  This  was  tone  by  either  interpolating 
or  extrapol  ati  ng,  in  time  and  space,  the  location  of  a  detected  gust  front 
such  that  it  coincided  with  the  location  of  a  particular  surface  site.  A  com- 
parvm  was  then  made  between  this  algorithm  indicated  time  of  passage  and  the 
a  c  1 1  a  1  time  of  passage  as  indicated  by  the  change  in  wind  direction  at  the 
surface  site.  The  results  are  as  follows: 

‘■lay  9,  19B1  case  (Figure  12) 

The  gust  front  algoritnm  was  run  for  three  times:  0006,  0013,  and 
0032  CST.  Unfortunately,  no  surface  stations  had  the  gust  front  pass  over 
them  from  0006  through  0032  CST.  Therefore,  in  order  to  do  a  comparison,  the 
position  of  the  gust  front  was  linearly  extrapolated  ahead  in  space  and  time 
to  the  nearest  three  stations.  The  results  are: 


Station 

number 


Actual  time  of  gust 
front  passage  (CST) 


Extrapolated  tine  of 
gust  front  passage  (CST) 


The  extrapolation  was  based  on  the  average  speed  of  the  gust  front  from  0006 
to  0032  CST.  The  wind  shift  at  all  three  stations  was  swift,  with  the  wind 
direction  going  from  SE  to  NW  in  one  minute  (the  lowest  time  scale  resolution 
for  NSSl  surface  stations  is  one  minute). 

April  26,  1934  case  (Figure  13) 

The  algorithm  was  run  for  three  times:  2037,  2042,  and  2047  CST.  For 
this  case,  two  surface  stations  are  located  along  the  detected  position  of  the 
gust  front  at  2037  CST.  The  time  of  passage  at  both  stations  was  around 
2035  CST.  The  time  of  passage  for  this  case  could  not  be  estimated  as  pre¬ 
cisely  as  in  the  previous  case  because  at  both  stations  it  took  about  four 
minutes  for  the  wind  direction  to  shift  from  S  to  WNW. 
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April  13,  *981  case  (Figure  14} 


Toe  clgoriton  was  again  run  for  three  times:  2111,  2116,  and  2123  CST. 
'nree  surface  stations  (No.  1-3)  are  located  on  jr  very  near  to  the  detected 
gust  front  positions.  In  this  case,  as  in  the  4/26/84  case,  the  surface  winds 
sniffed  gradually  from  a  SSE  direction  to  a  SW  direction,  and  finally  to  a  MW 


Fi;jur>:  14.  Same  as  Figure  1A ,  except  for  April  16,  1981. 


direction.  It  we  use  the  time  when  the  wind  direction  shifted  past  270°  for 
comparison,  we  get  time  differences  of  *1,  +1,  and  -2  minutes  for  the  three 
stations,  respectively,  with  a  positive  difference  meaning  the  algorithm  indi¬ 
cated  time  of  passage  was  ahead  of  the  surface  station  indicated  time  of 
passage,  -or  three  additional  sites  (No.  4-6),  we  can  extrapolate  ahead  in 
space  and  time  the  position  of  the  gust  front,  as  was  done  for  the  6/9/81 
case.  For  these  three  sites,  the  time  differences  are  -1,  -4,  and  -5  minutes, 
respecti vely . 

What  these  results  show  is  that,  in  most  cases,  there  was  fairly  good 
agreement  between  the  algorithm  output  and  the  surface  stations,  with  the  tin: 
difference  being  s  8  minutes  for  all  comparisons,  and  s  2  minutes  for  72%  of 
the  comparisons. 
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b.  Comparison  with  a  Second  Soppier  Radar 

A  second  comparison  of  algorithm  output  was  done  by  comparing  the  over¬ 
lapping  locations  of  gust  fronts  as  determined  from  two  different  Doppler 
radars  (the  NSSL  radars  located  at  Norman  and  Cimarron  (CIM),  wnich  is  about 
40  km  NW  of  Norman)  looking  at  the  same  gust  front.  The  comparison  was  made 
for  eight  different  elevation  scans  on  4/13/31,  with  the  gust  front  located 
between  10  and  30  km  away  from  each  radar  for  the  times  considered.  The 
results  are  given  in  Table  4.  (There  are  two  entries  at  2127  CST  and  0.4" 
becasue  of  a  gap  in  the  gust  fronts  detected  position,  as  seen  by  the  Nornar 
radar,  due  to  ground  clutter.)  The  root  mean  squared  (RMS)  distance  errors 
were  calculated  by  using  the  shortest  distance  between  any  two  points  of  the 
estimated  gust  front  locations.  We  see  that  the  average  RMS  distance  error  is 
less  than  .  km,  which  indicates  good  position  agreement.  Figure  15  snows  the 


algorithm  output  for  each  radar  at  2127  CST  for  approximately  the  same  ele¬ 
vation  angle.  From  Table  4  we  see  that  the  position  correlation  for  this  par¬ 
ticular  comparison  was  very  good,  with  an  RMS  distance  error  of  only  0.5  km. 


Figure  15.  Algorithm  output  for  two  different  radars;  a)  Norman  ('In 'O)  and 
o)  Cimarron  (CIM),  looking  at  the  name  gust  front  April  II, 
1981). 
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Results  of  the  Gust  Front  Position  Comparison 


Time  (CST) 

for  the  Norman 

Elevation  Angle  RMS 

and  Cimarron  Radars. 

distance  error  (km) 

No.  of  points 

2116 

O 

oc 

0 

1.4 

18 

2116 

1.3° 

1.1 

23 

2123 

0.4° 

0.7 

15 

2123 

0.9° 

0.6 

2/ 

2123 

1.3° 

0.7 

31 

2127 

0.4° 

1.4 

28 

2127 

0 

• 

o 

0.5 

12 

2127 

0.9° 

0.8 

76 

2127 

1.3° 

0.5 

72 

6.  SUMMARY 

Average 

=  0.88 

33.6 

The  gust  front  algorithm  originally  proposed  by  Uyeda  and  Zrnic'  (1985) 
was  fully  automated  and  made  ready  for  real-time  operational  testing.  Several 
changes  have  been  incorporated  in  order  to  simplify  the  algorithm,  and  speed 
up  its  running  time.  Of  substantial  impact  to  operational  users  is  our 
finding  that  it  may  be  necessary  to  operate  the  algorithm  with  more  than  one 
set  of  thresholds.  For  instance,  we  have  determined  that  higher  shear  (4  m 
s'1  km-*)  and  velocity  difference  (10  m  s'1)  thresholds  are  needed  for 
strongly-dynamic  weather  environments,  whereas  lower  thresholds  (2  m  s"1  km-1, 
5  m  s'1)  are  more  suited  for  dynamically  weaker  conditions.  All  tnree 
Oklahoma  gust  fronts  on  which  the  algorithm  was  tested  belonged  to  the  former 
category,  whereas  one  case  from  Denver  and  one  from  Memphis  were  associated 
with  weaker  weather  conditions.  When  appropriate  thresholds  were  used,  there 
were  very  few  false  alarms  (on  average  less  than  one  per  volume  scan).  But 
when  lower  thresholds  were  tried  for  one  of  the  Oklahoma  cases,  the  false 
alarm  rate  increased  substantially  (to  about  four  per  volume  scan).  Most  of 
these  were  due  to  problems  with  ground  clutter. 

Other  changes  to  the  algorithm  involved:  1)  the  inclusion  of  a  routine 
for  fitting  a  second-order  polynomial  to  the  gust  front  vector  positions  in 
cartesian  coordinates,  2)  inclusion  of  an  automated  vertical  continuity  check 


using  data  fr om  two  Tow-level  elevation  scans,  and  3)  inclusicn  of  a  sector¬ 
ized  uniform  wind  procedure  to  determine  winds  ahead  and  behind  the  front. 


The  outputs  of  the  algorithm  were  checked  for  consistency  and  compared  to 
data  from  other  sources.  For  instance,  vector  winds  from  the  uniform  wind  an¬ 
alysis  were  compared  with  surface  and  rawinsonde  observations,  as  well  as  with 
general  trends  of  the  Doppler  velocities.  Agreement  was  beyond  our  expec¬ 
tations;  moreover,  the  uniform  wind  analysis  proved,  on  at  least  one  occasion, 
to  be  a  better  indicator  of  gust  front  motion  than  tracking  of  the  centroid. 
Times  of  gust  front  passage  over  surface  stations  were  generally  within  two 
minutes  of  interpolated  or  extrapolated  positions  by  the  algorithm.  Also,  in 
the  only  i ntercompari son  between  the  positions  obtained  from  two  spaced 
Doppler  radars,  the  average  RMS  errors  were  less  than  D.9  km.  We  emphasize 
here  that  these  good  comparisons  are  for  strong,  single,  well-defined  gust 
fronts.  It  remains  to  be  determined  if  these  impressive  results  will  be  re¬ 
peated  in  situations  where  multiple,  weak  gust  fronts  are  prsent. 
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